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Abstract —  This  paper  considers  the  problem  of  detecting 
a  multichannel  signal  in  partially  homogeneous  environments, 
where  the  disturbances  in  both  test  signal  and  training  signals 
share  the  same  covariance  matrix  up  to  an  unknown  power  scal¬ 
ing  factor.  Two  different  parametric  Rao  tests,  referred  to  as  the 
normalized  parametric  Rao  (NPRao)  test  and  the  scale-invariant 
parametric  Rao  (SI-PRao)  test,  respectively,  are  developed  by 
modeling  the  disturbance  as  a  multichannel  autoregressive  (AR) 
process.  The  NPRao  and  SI-PRao  tests  entail  reduced  train¬ 
ing  requirement  and  computational  efficiency,  compared  with 
conventional  fully  adaptive,  covariance  matrix  based  solutions. 
The  SI-PRao  test  attains  asymptotically  a  constant  false  alarm 
rate  (CFAR)  that  is  independent  of  the  covariance  matrix  and 
power  scaling  factor  of  the  disturbance.  Comparisons  with  the 
covariance  matrix  based,  scale-invariant  generalized  likelihood 
ratio  test  (GLRT),  also  known  as  the  adaptive  coherent  estimator 
(ACE),  are  included.  Numerical  results  show  that  the  parametric 
Rao  detectors,  in  particular  the  SI-PRao  test,  attain  considerably 
better  detection  performance  and  use  significantly  less  training 
than  the  ACE  detector. 

Index  Terms — Rao  test,  parametric  models,  partially  homo¬ 
geneous  environment,  multichannel  signal  detection,  space-time 
adaptive  signal  processing. 

I.  Introduction 

Multichannel  adaptive  detection  is  encountered  in  numerous 
applications  including  radar  [1],  [2],  wireless  communications 
[3],  hyperspectral  imaging  [4],  [5],  and  others.  Space-time 
adaptive  processing  (STAP)-based  multichannel  adaptive  de¬ 
tection  has  been  successfully  utilized  to  mitigate  the  effect 
of  clutter  and/or  interference  in  radar,  remote  sensing,  and 
communication  systems  [1] — [5]. 

Traditional  STAP  detectors  are  developed  usually  for  homo¬ 
geneous  environments,  where  the  disturbances  in  both  the  test 
and  training  signals  are  assumed  to  be  independent  and  iden¬ 
tically  distributed  (i.i.d.).  Examples  include  the  Reed,  Mallet, 
and  Brennan  detector  [6],  Kelly’s  generalized  likelihood  ratio 
test  (GLRT)  [7],  the  adaptive  matched  filter  (AMF)  detector 
[8],  [9],  Rao  test  [10],  among  others.  For  the  above  STAP 
detectors,  estimation  of  the  space-time  covariance  matrix  from 
training  signals  requires  a  large  number  of  training  signals 
and  excessive  computation  power,  especially  when  the  joint 
space-time  dimension  is  large.  To  alleviate  these  problems, 
parametric  STAP  detectors  have  been  developed  by  modeling 
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the  disturbance  as  a  multichannel  autoregressive  (AR)  process, 
such  as  the  parametric  AMF  (PAMF)  [11],  parametric  Rao  test 
[12],  and  parametric  GLRT  [13]. 

In  this  paper,  we  consider  multichannel  signal  detection 
in  partially  homogeneous  environments,  where  the  test  signal 
shares  the  same  covariance  matrix  with  the  training  signals  up 
to  an  unknown  power  scaling  factor  under  the  null  hypoth¬ 
esis  [  14]— [  17] .  This  scenario  is  motivated  by  the  following 
observation:  A  number  of  guard  cells  used  in  the  STAP 
implementation  to  separate  the  test  signal  and  training  signals 
may  lead  to  a  power  difference  between  the  test  and  training 
signals.  Specifically,  the  training  signals  over  range  cells  are 
assumed  to  be  i.i.d.  CM( 0,  R)  which  denotes  a  complex 
Gaussian  distribution  with  zero  mean  and  covariance  matrix 
R,  while  the  disturbance  in  the  test  signal  is  independent 
with  the  training  signals  with  distribution  CM (0,  AR),  where 
A  denotes  the  unknown  power  scaling  factor.  Depending  on 
a  priori  knowledge  about  the  power  scaling  factor  A,  or 
the  covariance  matrix  R,  or  both,  non-adaptive  and  adap¬ 
tive  STAP  detectors  have  been  developed  in  [  15]— [  1 8].  The 
scale-invariant  GLRT,  which  is  also  known  as  the  adaptive 
coherence  estimator  (ACE),  is  first  introduced  in  [18]  for  the 
compound  Gaussian  environment  and  also  developed  in  [16] 
for  the  partially  homogeneous  environment.  It  is  shown  that 
the  ACE  coincides  with  the  Rao  and  Wald  tests  in  the  partially 
homogeneous  environment  [17],  which  is  also  a  constant  false 
alarm  rate  (CFAR)  detector.  As  a  covariance  matrix  based 
STAP  detector,  the  ACE  needs  to  estimate  the  covariance 
matrix  from  the  target-free  training  signals  which  entails  large 
training  requirement,  and  to  invert  the  covariance  matrix  for 
implementation  also  leads  to  a  high  computational  complexity. 
At  a  minimum,  K  >  JN  i.i.d.  training  signals  are  needed  to 
ensure  a  full-rank  estimate  of  the  JN  x  JN  covariance  matrix 
R,  where  J  denotes  the  number  of  spatial  channels  and  N  the 
number  of  temporal  observations.  For  example,  the  KASSPER 
dataset  consists  of  J  =  11  spatial  channels  and  N  =  32 
coherent  pulses  for  a  total  of  JN  =  352  spatial-temporal 
dimension.  Such  a  demanding  training  requirement  usually 
cannot  be  met  in  practice.  It  also  makes  the  homogeneous 
(i.e.,  i.i.d.)  assumption  across  such  a  broad  range  of  training 
cells  impractical. 

To  address  this  issue,  we  take  an  approach  where  the  distur¬ 
bance  in  the  partially  homogeneous  environment  is  modeled 
by  a  multichannel  AR  process.  This  modification  results  in 
a  parametric  STAP  detector  which  requires  a  relaxed  local 
homogeneous  assumption  for  the  training  signals:  that  is,  the 
training  signals  are  assumed  to  be  i.i.d.  over  a  small  duration, 
as  opposed  to  the  excessive  homogeneous  training  signals 
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required  by  the  ACE.  The  local  homogeneous  assumption  is 
reasonable  because  the  effect  of  clutter  variation  across  a  small 
area  is  generally  negligible.  In  the  compound  Gaussian  model 
[18]— [23],  it  is  often  assumed  that  the  texture  component 
is  slow  varying,  which  is  similar  to  the  local  homogeneous 
assumption. 

In  this  paper,  parametric  detection  by  using  the  Rao  prin¬ 
ciple  is  considered  for  the  partially  homogeneous  environ¬ 
ment.  The  first  parametric  Rao  test,  which  is  referred  to  the 
normalized  parametric  Rao  (NPRao)  test,  is  developed  in  a 
heuristic  way.  Specifically,  it  first  assumes  the  knowledge  of 
the  power  scaling  factor,  applies  the  Rao  principle  and  obtains 
a  Rao  test  that  depends  on  the  power  scaling  factor,  and  finally 
replaces  the  unknown  power  scaling  factor  by  some  estimate. 
The  second  Rao  test,  referred  to  the  scale-invariant  parametric 
Rao  (SI-PRao)  test,  is  developed  by  applying  the  Rao  principle 
with  the  joint  maximum  likelihood  (ML)  estimates  of  all 
unknown  parameters  in  the  null  hypothesis.  The  asymptotical 
distribution  of  the  SI-PRao  test  is  derived  in  closed  form.  It  is 
shown  that  the  SI-PRao  test  is  asymptotically  independent  of 
the  unknown  parameters  in  the  null  hypothesis,  which  results 
in  the  property  of  CFAR.  Comparisons  between  the  ACE 
detector  show  that  our  parametric  Rao  tests  in  the  partially 
homogeneous  environment  exhibit  improved  detection  perfor¬ 
mance  when  the  training  signals  are  limited,  i.e.,  when  K 
is  too  small  to  yield  a  reliable  estimate  of  the  space-time 
covariance  matrix  R  that  is  required  by  the  ACE  detector. 

The  rest  of  this  paper  is  organized  as  follows.  Section  II 
contains  the  data  model  and  problem  statement.  Prior  solutions 
to  the  problem  of  interest  are  briefly  reviewed  in  Section  III. 
The  NPRao  and  SI-PRao  tests  are  developed  and  summarized 
in  Section  IV.  Asymptotical  performance  of  the  SI-PRao 
detector  is  also  included  in  this  section.  Numerical  results 
are  presented  in  Section  V.  Finally,  Section  VI  contains  the 
concluding  remarks. 

II.  Data  Model  and  Problem  Statement 
A.  Data  Model 

Consider  the  problem  of  detecting  a  known  multichannel 
signal  with  unknown  amplitude  in  the  presence  of  spatially 
and  temporally  correlated  disturbance:  (e.g.,  [1]): 

H0-  x0(n)  =  d0(n),n  =  0,1,-  •  •  ,N  -  1 

Hi  :  xo (n)  =  as(n)  +  do (n),n  =  0, 1,  -  -  -  ,  N  —  1 

where  all  vectors  are  J  x  1  vectors,  J  denotes  the  number 
of  spatial  channels,  and  N  is  the  number  of  temporal  obser¬ 
vations.  In  the  sequel,  xo(n)  is  referred  to  as  the  test  signal, 
s(n)  is  the  signal  to  be  detected  with  amplitude  a,  and  do(n) 
is  the  disturbance  signal  that  may  be  correlated  in  space  and 
time.  Besides  the  test  signal  xo(n),  there  may  be  a  set  of 
training  signals  Xfc(n),  k  =  1, 2,  •  •  •  ,  K,  to  assist  in  the  signal 
detection: 

Xfc(n)  =  dfc(n),  n  =  0, 1,  •  •  •  ,  V  —  1.  (2) 

In  radar  systems,  training  signals  may  be  obtained  from  range 
cells  adjacent  to  the  test  cell.  However,  training  signal  is 


generally  limited  or  may  even  be  unavailable.  In  the  extreme 
training-free  case,  we  have  K  =  0. 

Define  the  following  JN  x  1  space-time  vectors: 

s  =  [sT(0),sT(1),  •  •  •  ,sT(N  —  1)]T  , 

dfc  =  [dfc  (0),  d^(l),  •  •  •  ,dfr(N  —  1)]T  , 

xfc  =  [x^  (0),  xjT (1),  •  •  •  ,xl(N-  1)]T, 

where  k  =  0, 1,  •  •  •  ,  K.  It  follows  that  (1)  can  be  rewritten  in 
a  compact  form 

Hq  :  xo  =  do 

(3) 

Hi  :  x0  =  as  +  d0. 

The  binary  composite  hypothesis  testing  problem  is  to  select 
between  Hq  :  a  =  0  and  Hi  :  a  ^  0. 

The  general  assumptions  for  the  STAP  in  partially  homo¬ 
geneous  environment  are  [  15]— [  17],  [24],  [25]: 

•  AS1:  The  signal  vector  s  is  deterministic  and  known  to 
the  detector; 

•  AS2:  The  signal  amplitude  a  is  complex-valued,  deter¬ 
ministic,  and  unknown; 

•  AS3:  The  disturbance  signals  do  and  {&k}k=i  are 
mutually  independent  with  distribution  CM ( 0 ,  AR)  and 
CM (0,  R),  respectively,  where  A  >  0  is  an  unknown 
power  scaling  factor. 

If  A  =  1,  the  partially  homogeneous  environment  reduces 
to  the  homogeneous  case.  In  this  paper,  a  multichannel  AR 
process  is  employed  to  model  the  disturbance  and  the  AS3  is 
modified  as  follows 

•  AS3':  The  disturbance  signal  dfc(n),  k  =  0, . .  .,K,  can 
be  modeled  as  a  J-channel  AR(P)  process  with  model 
order  P: 

dfc(n)  =  -  Yl?=i  Aff  (i)dfe(n  -  i)  +  £fc(n),  (4) 

where  {Ai/(i)}^1  denote  the  unknown  J  x  J  AR 
coefficient  matrices,  £k(ji)  denote  the  J  x  1  spatial  noise 
vectors  that  are  temporally  white  but  spatially  colored: 
e0(n)  ~  CM(0,  AQ)  and  {efe(n)}^=1  ~  CN(0,  Q), 
respectively,  where  A  is  the  unknown  power  scaling  factor, 
and  Q  denotes  the  unknown  J  x  J  spatial  covariance 
matrix. 

Note  that  the  power  scaling  factor  A  on  the  spatial  covariance 
matrix  ensures  the  same  power  scaling  factor  on  the  spatial- 
temporal  covariance  matrix.  The  problem  of  interest  here  is 
to  develop  a  decision  rule  for  the  above  composite  hypothesis 
testing  problem  based  on  assumptions  AS1,  AS2,  and  AS3'. 

III.  Prior  Solutions 

Depending  on  the  amount  of  a  prior  knowledge  about  the 
unknown  parameters,  a  number  of  solutions  for  the  STAP 
detection  in  partially  homogeneous  environments  have  been 
proposed.  If  the  space-time  covariance  matrix  R  and  the  power 
scaling  factor  A  are  both  known  exactly,  the  optimal  detector 
is  the  phase  invariant  matched  filter  (PIMF)  [20],  [25] 

|s-ffR_1X0|2  H1 

Tpimf  =  — — 77— — - —  ^  7pimf,  (5) 

As  R_is  H0 
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where  7pimf  denotes  the  PIMF  threshold  subject  to  a  selected 
probability  of  false  alarm.  In  the  case  of  unknown  A,  the 
normalized  match  filter  (NMF)  replaces  A  with  an  estimate, 
which  is  given  by 


Tnmf  — 


|shR  1Xo| 


(s^R^s)  (x^R^xq)  h0 


Ho 

2;  7nmf, 


(6) 


where  7nmf  denotes  the  NMF  threshold.  It  should  be  noted 
that  both  the  non-adaptive  PIMF  and  NMF  detectors  cannot 
be  implemented  in  practice  since  R  is  unknown,  but  they  can 
be  used  as  a  baseline  for  performance  comparison. 

Adaptive  detectors  are  formulated  by  estimating  the  covari¬ 
ance  matrix  R  from  the  target-free  training  signals 


1  K 


xfcxfc  ■ 


(7) 


fc= l 


The  normalized  adaptive  matched  filter  (NAMF)  [18],  also 
known  as  the  adaptive  coherence  estimator  (ACE)  [15],  is 


Pace 


sffR 

lx0 

2 

1  1 

^x^R-lxo) 

Ho 

2;  7ace, 
Ho 


(8) 


where  7ace  denotes  the  ACE  threshold.  It  is  shown  that  the 
ACE  is  equivalent  to  the  scale-invariant  GLRT  in  the  partially 
homogeneous  environment  [15]. 

By  utilizing  a  multichannel  AR  model,  the  normalized 
parametric  adaptive  matched  filter  (NPAMF)  is  reported  [20], 
[21],  [26]: 


A.  Normalized  Parametric  Rao  Test 

The  NPRao  test  is  a  heuristic  detector  obtained  via  a  two- 
step  approach.  First,  by  assuming  that  the  power  scaling 
parameter  A  is  known,  a  Rao  test  can  be  obtained  by  finding 
the  ML  estimates  of  the  AR  coefficient  matrices  {A H  (p)} 
and  spatial  covariance  matrix  Q.  The  derivation  is  similar  to 
the  Rao  test  in  [12]  for  a  homogeneous  environment.  Hence, 
we  skip  the  details  and  just  present  the  final  result.  The  new 
detector,  conditioned  on  a  known  A,  is  given  by 


where  s  and  xp  are  the  temporally  whitened  steering  vector 
and  test  signal,  respectively,  by  using  a  set  of  estimates 
{AH(p)}  of  the  AR  coefficient  matrices: 

p 

s  (n)  =s  (n)  +  AH  {p)s  (n  -  p) ,  (11) 

P=  1 
P 

x0  (n)  =xo  (n)  +  ^2AH  (p)x0  (n-p),  (12) 

p=  i 


where  AH 

[12]: 


Aff(l),Aff(2),---  ,Ah  (P) 


is  given  by 


A  -H  _  1 

A  ^yx^yy  ' 


(13) 


Tnpamf  — 

n~1~h  „  „  2 

E  §p  («)Qp  : 1*o,P  (n) 

n=P 

N—  1  A H  ~  ~  \  /TV- 1  I  I 

E  §p  (n)Qp1sP(?r)  E  xo,p(n)Qp  1xo,p  («) 
n=P  J  \n—P  J 

Hx 

2;  Tnpamf,  (9) 

Ho 

where  7npamf  denotes  the  NPAMF  threshold,  Qp  denotes  an 
estimate  of  the  spatial  covariance  matrix  Q,  xp.p  (n)  and  §p  (n) 
are  the  temporally  whitened  test  and  steering  vector,  respec¬ 
tively.  The  NPAMF  is  originally  developed  as  a  solution  for 
detection  in  the  compound-Gaussian  environment  [20],  [21]; 
it  can  also  be  used  in  the  partially  homogeneous  environment. 


The  Q  in  (10)  is  an  estimate  of  the  spatial  covariance  matrix 

Q  [12]: 

Q  =  ^  ~  >  (14) 

where  L  =  ( K  +  1  ){N  —  P)  and 

K  N-l 

Rxx  =^2  Xfc  (n)x^  (n)>  (15) 

k= 0 n—P 
K  N-l 

&yy  yk  (n) ,  (16) 

k—0 n—P 
K  N-l 

Ryz  =  yk  (n)  {n),  (17) 

k—0 n—P 


IV.  Parametric  Rao  Test 
In  Partially  homogeneous  Environment 

In  this  section,  two  parametric  Rao  detectors  based  on 
assumptions  AS1,  AS2,  and  AS3'  are  developed.  The  NPRao 
test  is  developed  in  a  manner  similar  to  the  NPAMF,  while 
the  SI-PRao  test  is  derived  by  finding  the  ML  estimates  of  all 
nuisance  parameters  under  the  null  hypothesis. 


and  yjt(n)  =  [x£(n  —  1),  ...,  x^(n- P)]T  G  CJPxl1k  = 

0 

Second,  we  replace  the  power  scaling  factor  A  in  (10)  by 
some  estimate.  One  candidate  is  given  by  [20],  [21]  (also  see 
(9)) 

IV-l 

A  =  ^  xf  (n)Q_1x0  (n).  (18) 

n—P 
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This  leads  to  our  NPRao  detector  given  by 


where  7NPRao  denotes  the  NPRao  test  threshold. 

It  is  noted  that  the  NPAMF  detector  is  equivalent  to  the 
NPRao  test  with  one  exception:  while  the  NPAMF  detector 
uses  only  training  signals  for  parameter  estimation,  the  NPRao 
detector  uses  both  training  and  test  signals  for  estimation.  The 
difference  is  similar  to  the  one  discussed  in  [12]  between  the 
PAMF  detector  and  the  parametric  Rao  detector. 


J=4 


Fig.  1.  Numerical  evaluation  of  /(A)  for  several  sets  of  system  parameters 
when  J  =  4  and  A  =  4 


~  _!  ~  -  _i  ~  ~  -  _I 

matrices  R0  2  R/yR0  2  and  R0  2  RjfiSR0  2 ,  respectively,  and 


B.  Scale-Invariant  Parametric  Rao  Test 


The  heuristic  NPRao  test  cannot  ensure  the  invariance  to 
the  power  scaling  factor,  which  is  needed  for  effective  clutter 
mitigation  in  the  partially  homogeneous  environment  [15], 
[16].  To  this  end,  a  parametric  Rao  test  by  evaluating  the  ML 
estimates  of  all  nuisance  parameters  including  A,  A  and  Q 
is  developed  and  the  resulting  detector,  the  SI-PRao  test,  is 
shown  to  be  invariant  to  the  power  scaling  factor. 

Specifically,  as  shown  in  Appendix,  the  SI-PRao  test  is 
given  by 


^SI-PRao 


NE/(n-rA)Q-i(x)*0(n-A) 


X^/{n-,X)Q-^x)^nrx) 


s;7si-PRao,  (20) 

H0 


The  A  in  (20)  is  the  ML  estimate  of  A,  which  is  obtained  by 
solving  the  following  equation: 


J 

K  +  1 


J(p+ i)  1  Jp  i 

1  +  A  ^  ^  1  +  Xvi 

1=1  1=1 


(21) 


where  f+^  and  are  the  eigenvalues  of  the 


Rn  — 


N-l  N- 1 

E  yo(n)y^(n)  E  yo  (n)x^(n) 

n=P  n=P 

N-l  N-l 

E  x0(n)y^(ra)  E  x0(n)x^(n) 

_  n=P  n=P 


(22) 


Ra'  = 


K  N-l 


K  N-l 


EEwWyfW  EEyfcWxf(n) 


k= 1 n=P 
K  N-l 


k= 1 n=P 
K  N-l 


E  E  xfc  (n)  yf  (n)  EExfc(n)xf(n) 


k=l n=P 


k=l n=P 


N-l 


Ro,y  =  yo(n)yo  (n)> 

n=P 
K  N-l 

R/r,y  =  yk  (n)- 


(23) 

(24) 

(25) 


k= 1 n=P 


Remark:  The  existence  of  a  non-negative  solution  A  to 

_ 1_  ^  ^  _1 

(21)  is  guaranteed.  To  see  this,  note  that  R0  2  R/yR0  2  and 
R-0,y  are  positive  definite.  Therefore,  /i,;  and  are 

positive.  The  function 


/(A) 


J 

K+  1 


J(P+ 1)  1  Jp  i 

^  1  +  A  Hi  +  ^  1  +  Xvi 
2=1  2=1 


(26) 


is  a  continuous  function  on  (0, +oo).  When  K  >  0,  i.e.,  at 
least  one  training  signal  is  available,  we  have 

lim /(A)  =  — -  J  <  0,  (27) 

lim  /(A)  =  J  >  0,  (28) 

A— >oo  A  +  1 

which  implies  that  there  is  at  least  one  A  within  in  the  interval 
(0,  Too)  giving  /(A)  =  0.  On  the  other  hand,  the  uniqueness 
of  the  ML  estimate  of  A  is  more  difficult  to  establish  due 
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to  the  nonlinear  dependence  on  the  parameters  p.-t  and  //, . 
Numerical  examples  suggest  that  the  solution  is  unique.  For 
illustration  purpose.  Figure  1  shows  the  numerical  evaluation 
of  /(A)  for  several  sets  of  system  parameters.  It  is  seen  that 
the  function  /(A)  is  a  monotonically  increasing  function  over 
A  €  (0,+oo),  which,  combining  with  (27)  and  (28),  suggests 
that  the  solution  is  unique.  However,  a  more  rigorous  proof  is 
not  available. 

Once  the  ML  estimate  of  A,  the  temporally  whitened  steer¬ 
ing  vector  s [n\  A)  and  test  signal  xo(n;  A)  in  (20)  are  obtained 
from  (11)  and  (12)  with  the  following  ML  estimate  of  A 

Ah  =  -R"  (A)  R-1  (a)  .  (29) 

Finally,  Q(A)  in  (20),  the  ML  estimate  of  Q,  is  given  by 

Q  (a)  =  \  (Rxx(A)  -  R"(A)R,-;(A)R?/,(A))  ,  (30) 

where 

TV— 1  K  TV— 1 

Rxx(A)  =  A-1  Y,  xo  ( n )  x^  (n)  +  Y  Y  Xfc  (n)  Xk  (n)’ 

n—P  k—1 n—P 

(31) 

N-l  K  N-l 

Ry?/(A)  =  A-1  Y  yo  (n)  Yo  (n)  +  YYl  yk  (n)  yk  (”). 

n—P  k—1 n—P 

(32) 

JV-l  K  N-l 

&yx(  A)  =  A'1  Y  yo  ( n )  *0  H  +  YY1  yk  M  («)• 

n—P  k—1 n—P 

(33) 


C.  Asymptotic  Detection  Performance 


According  to  [27],  the  asymptotic  distribution  of  the  Rao 
test  statistics  can  be  obtained  as 


a  (  xl  under  H0 
SI'PRao  l  X?  (P)  ■  under  Hi 


(34) 


where  x‘2  denotes  the  central  Chi-squared  distribution  with  2 
degrees  of  freedom  and  X2  (p)  the  noncentral  Chi-squared 
distribution  with  2  degrees  of  freedom  and  noncentrality 
parameter  p 


n  N~1 

2|a|2  sff(n)Q^1s(n) 

n—P 


where  s  is  the  temporally  whitened  steering  vector  given  by 
(11)  but  with  A  replaced  by  the  true  AR  coefficient  matrix 
A.  From  (35),  it  is  ready  to  show  that,  for  a  given  threshold, 
the  probability  of  false  alarm  is 

D  (  7SI-PRao\ 

Pf  =  exp  - - — j  ,  (36) 

which  shows  that  the  statistic  of  the  SI-PRao  test  under  //0  is 
independent  of  the  power  scaling  factor  and  the  covariance 
matrix,  and  further  implies  the  SI-PRao  test  is  a  CFAR 
detector,  while  the  probability  of  detection  is 

Pd=  [  ^exp  I0{Vpx)dx,  (37) 

‘'TSI-PRao  ^  \  ^  / 

where  Iq  (x)  is  the  modified  Bessel  function  of  the  first  kind 
and  zeroth  order. 


K=2,  J=4,  N=32,  P=0.01,  AR(2),  A=4 


Fig.  2.  Probability  of  detection  versus  SINR  when  P&  =  0.01,  J  =  4, 
TV  =  32,  K  =  2,  and  A  =  4. 


V.  Numerical  Examples 

We  now  report  simulation  results  for  the  proposed  detectors. 
Throughout  this  section,  the  disturbance  signal  is  generated 
as  a  multichannel  AR(2)  process  with  AR  coefficient  A  and 
a  spatial  covariance  matrix  Q.  These  parameters  are  set  to 
ensure  that  the  AR  process  is  stable  and  Q  is  a  valid  covariance 
matrix,  but  otherwise  are  randomly  selected.  The  signal  vector 
s  corresponds  to  a  uniform  equispaced  linear  array  with  ,7  =  4 
antenna  elements,  N  temporal  pulses,  and  randomly  selected 
normalized  spatial  frequency  u>s  and  Doppler  frequency  tod- 
The  steering  vector  is  given  by 

s  =  st(wd)  (gis s(u>s),  (38) 

where  s t(u>d)  denotes  the  N  x  1  temporal  steering  vector 

stM  =  -±=[l,e^V--  (39) 

and  ss(ws)  denotes  the  J  x  1  spatial  steering  vector 

ss(Ws)  =  -i=M"v  (40) 

The  signal-to-interference-plus-noise  ratio  (SINR)  is  defined 
as 

SINR  =  |a|2sffR-1s,  (41) 

where  the  JN  x  JN  covariance  matrix  R  can  be  uniquely 
determined  once  A  and  Q  are  selected. 

Numerical  simulations  are  provided  for  various  values  of 
K,  N,  and  A  when  ,7  =  4  and  P  =  2.  In  particular,  we 
consider  two  distinct  cases:  1)  the  limited-training  case,  e.g., 
K  =  2  and  I\  =  4;  and  2)  the  asymptotic  case,  e.g.,  K  =  64. 
The  simulation  results  are  shown  in  terms  of  the  probability  of 
detection  versus  SINR  when  the  probability  of  false  alarm  is 
fixed  as  Pf  =  0.01.  In  each  case,  the  optimum  but  non-adaptive 
PIMF  detector  (5),  whose  Pd  and  Pf  can  be  determined 
analytically  [16],  is  included  as  a  benchmark  (note  that  the 
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K=2,  J=4,  N=32,  P=0.01,  AR(2),  A=8 


0.9 


0.8 


:  0.7  - 


S  0.6  - 


;  0.5  - 


i0-4 

CO 

o  0.3 

£ 

0.2- 


0.1  - 


- PIMF  (analytic) 

O  SI-PRao  (asymptotic) 
+  SI-PRao  (simulation) 
x  NPRao  (simulation) 
ACE  (K=2JN,  analytic) 


- r[y(^  0  0  4)  O 

P+  '  x 

/■+  ;-x  x 


<6 

0  +  , 

0  .  x. 


0  _ 


/  + 
0 


0 

/+ 

P  X  / 

'  x  / 

p+  , 

+  X 


5  10  15 

SINR  (dB) 


20 


25 


K=4,  J=4,  N=16,  P=0.01,  AR(2),  A=8 


Fig.  3.  Probability  of  detection  versus  SINR  when  Pi  =  0.01,  J  =  4,  Fig.  5.  Probability  of  detection  versus  SINR  when  P\  =  0.01,  J  =  4, 
N  =  32,  K  =  2,  and  A  =  8.  N  =  16,  K  =  4,  and  A  =  8. 


K=4,  J=4,  N=16,  P  =0.01,  AR(2),  A=4 


Fig.  4.  Probability  of  detection  versus  SINR  when  Pi  =  0.01,  J  =  4, 
At  =  16,  A"  =  4,  and  A  =  4. 


PIMF  is  a  clairvoyant  detector  that  has  full  knowledge  of  all 
unknown  parameters).  The  ACE  detector,  whose  performance 
can  also  be  determined  analytically  (see  [16],  [28]),  is  also 
included  to  show  how  the  parametric  detectors  compare  with 
this  popular  covariance  matrix  based  detector  in  partially  non- 
homogeneous  environments.  Since  the  ACE  detector  needs  a 
full-rank  space-time  covariance  matrix  estimate,  we  assumes  it 
has  sufficiently  i.i.d.  training  signals  with  K  =  2  JN.  Finally, 
we  include  the  simulated  results  for  the  proposed  NPRao  and 
SI-PRao  detectors;  for  the  latter,  we  also  include  its  asymptotic 
performance  as  shown  in  Section  IV-C. 


K=2,  J=4,  N=32,  AR(2),  A=4 
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Fig.  6.  Comparison  of  SINR  loss  of  the  SI-PRao  detector  with  respect  to  the 
PIMF  in  cases  of  Pf  =  (0.1,  0.01,  0.001}when  J  =  4,  N  =  32,  K  =  2, 
P  =  2,  and  A  =  4. 


A.  Limited-Training  Case 

In  this  case,  we  first  consider  two  scenarios:  1)  K  =  2 
and  N  =  32,  where  the  training  signals  are  limited  and  the 
number  of  temporal  samples  is  moderate;  2)  K  =  4  and 
N  =  16,  where  we  have  slightly  more  training  signals  but 
the  temporal  samples  are  limited.  For  K  =  2  and  N  =  32, 
Figs.  2  and  3  show  the  probability  of  detection  versus  SINR 
when  A  =  4  and  A  =  8,  respectively.  It  is  seen  that,  for  a  small 
value  of  A  =  4,  the  NPRao  and  SI-PRao  tests  with  K  =  2 
generally  outperform  the  ACE  detector  with  K  =  256.  For  a 
larger  value  of  A  =  8  as  shown  in  Fig.  3,  the  SI-PRao  test 
behaves  similarly  to  the  PIMF,  while  the  NPRao  test  degrades 
considerably,  indicating  that  it  is  not  invariant  to  the  power 


7 


K=2,  J=4,  N=32,  AR(2),  A=4  K=64,  J=4,  N=1 6,  P  =0.01 ,  AR(2),  A=8 


SINR(dB)  SINR(dB) 


Fig.  7.  Comparison  of  SINR  loss  of  the  NPRao  detector  with  respect  to  the  Fig.  9.  Probability  of  detection  versus  SINR  when  PA  =  0.01,  J  =  4, 

PIMF  in  cases  of  Pf  =  {0.1,  0.01,  0.001}when  J  =  4,  N  =  32,  K  =  2,  TV  =  16,  K  =  64,  and  A  =  8. 

P  =  2,  and  A  =  4. 


SINR  (dB) 


Fig.  10.  Effect  of  the  number  of  pulses  on  the  detection  performance  of  the 
Fig.  8.  Probability  of  detection  versus  SINR  when  Pi  =  0.01,  J  =  4,  SI-PRao  detector  when  Pf  =  0.01.  J  =  4,  K  =  2,  P  =  2  and  A  =  4. 

N  =  16,  if  =  64,  and  A  =  4. 


scaling  factor  A. 

In  the  second  case  with  K  =  4  and  N  =  16,  the  simulation 
results  are  plotted  in  Figs.  4  and  5,  respectively,  for  A  =  4 
and  A  =  8.  It  is  seen  that,  the  SI-PRao  test  yields  the  best 
performance  among  the  adaptive  detectors.  Both  the  SI-PRao 
test  and  ACE  detector  are  seen  to  be  invariant  to  the  power 
scaling  factor.  However,  the  NPRao  test  shows  much  worse 
detection  performance  than  the  SI-PRao  test  and  the  ACE 
detectors  for  both  A  =  4  and  A  =  8. 

We  next  examine  the  SINR  loss  of  the  proposed  Rao 
detectors  with  respect  to  the  PIMF  in  cases  of  different  proba¬ 
bilities  of  false  alarm.  The  results  of  the  SI-PRao  and  NPRao 
detectors  are  shown  in  Fig.  6  and  Fig.  7,  respectively,  when 
several  probabilities  of  false  alarm  Pf  =  {0.1,0.01,0.001} 


are  considered.  It  is  noticed  that,  for  a  given  probability  of 
detection,  both  Rao  detectors  show  a  larger  SINR  loss  with 
respect  to  the  PIMF  when  the  probability  of  false  alarm  is 
smaller.  However,  the  performance  loss  of  the  SI-PRao  test  is 
generally  smaller  than  the  NPRao  test  in  all  three  cases. 

B.  Asymptotic  Case 

An  asymptotic  scenario  with  large  I\  is  simulated  to  verify 
the  asymptotic  performance  of  the  SI-PRao  test  derived  in 
Section  IV-C.  The  simulation  parameters  are  J  =  4,  N  =  16, 
P  =  2,  and  K  =  JN  =  64.  The  results  are  shown  in  Figs.  8 
and  9  with  A  =  4  and  A  =  8,  respectively.  It  is  seen  that, 
the  probability  of  detection  obtained  by  simulation  approaches 
the  asymptotic  performance  of  the  SI-PRao  test  in  both  cases. 
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Fig.  11.  Probability  of  detection  versus  SINR  when  model  order  of 
multichannel  AR  process  used  for  computing  test  statistic  is  true  (P  =  2), 
overestimated  (assuming  P  =  3),  and  underestimated  (assuming  P  =  1), 
when  Pf  =  0.01,  J  =  4,  N  =  32,  K  =  8,  and  A  =  4. 

Also,  with  a  plenty  of  training  signals,  the  NPRao  test  provides 
detection  performance  almost  identical  to  that  of  the  SI-PRao 
test.  The  performance  of  the  ACE  detector  with  K  =  128, 
which  is  twice  as  that  of  the  Rao  tests,  is  also  included  for 
comparison. 

It  is  also  interesting  to  show  the  impact  of  the  sample 
size,  i.e.,  the  number  of  pulses  number  N,  on  the  detection 
performance.  The  results  are  shown  in  Fig.  10,  where  the 
detection  probability  of  the  SI-PRao  detector  converges  as 
N  increases  and  approaches  the  asymptotic  performance.  The 
NPRao  detector  has  a  similar  behavior,  which  is  not  shown 
for  space  limitation. 

C.  Effect  of  Model  Mismatch 

The  above  simulation  examples  are  based  on  two  assump¬ 
tions:  1)  the  model  order  is  known;  2)  the  disturbance  is 
exactly  a  multichannel  AR  process.  In  this  section,  we  evaluate 
the  detection  performance  of  the  proposed  Rao  detectors  when 
these  assumptions  are  not  met. 

We  first  consider  the  case  when  the  disturbance  is  an  AR 
process,  but  there  is  an  model  estimation  error.  In  practice, 
one  needs  to  estimate  the  model  order  P  of  the  multichannel 
AR  before  application  of  the  SI-PRao  and  NPRao  detectors.  A 
model  order  estimation  procedure  may  yield  a  small  estimation 
error.  Fig.  1 1  depicts  the  detection  performance  of  the  SI-PRao 
detector  when  the  model  order  is  under-  and  over-estimated, 
respectively.  It  is  seen  that  the  order  mismatch  causes  some 
performance  degradation.  However,  in  both  cases,  the  degrada¬ 
tion  is  not  significant.  It  is  also  seen  that  overestimation  of  the 
model  order  has  a  smaller  effect  on  the  detection  performance 
than  underestimation.  This  behavior  is  similar  to  the  case  of 
the  standard  parametric  Rao  test  in  [12]. 

Next,  we  consider  the  case  when  the  disturbance  is  not  a 
multichannel  AR  process.  To  show  the  detection  performance 


Fig.  12.  Probability  of  detection  versus  SINR  for  for  the  KASSPER  dataset 
when  Pf  =  0.01,  J  =  11,  N  =  32,  K  =  8,  P  as  1  and  A  =  4. 

in  a  more  realistic  environment,  we  use  the  KASSPER  dataset, 
which  contains  many  challenging  real-world  effects,  including 
heterogeneous  terrain,  array  errors,  and  dense  ground  targets 
(see  [29]  for  a  detailed  description  of  the  KASSPER  dataset). 

Fig.  12  shows  the  probability  of  detection  versus  SINR  for 
the  KASSPER  dataset  when  Pf  =  0.01,  J  =  11,  N  =  32, 
K  =  8,  P  =  1  and  A  =  4.  The  covariance  matrix  of  the  test 
signal  corresponds  to  range  cell  r200,  whereas  the  covariance 
matrices  of  two  ( K  =  2)  training  signals  are  from  range 
cells  rl97  and  r203,  i.e.,  with  two  guard  cells  between  the 
test  cell  and  training  cells.  This  scenario  includes  some  non- 
homogeneous  effect  since  the  covariance  matrices  for  the  test 
and  training  cells  are  different.  As  shown  in  Fig.  12,  with 
two  training  signals  and  by  modeling  the  disturbance  as  an 
AR(1)  process,  the  SI-PRao  and  NPRao  detectors  achieve  a 
close  performance  that  is  only  2  to  3  dB  away  from  the  PIMF 
detector. 

VI.  Conclusion 

The  multichannel  adaptive  detection  in  the  partially  ho¬ 
mogeneous  environment  by  modeling  the  disturbance  as  a 
multichannel  AR  process  has  been  addressed.  Two  parametric 
Rao  tests  have  been  developed  to  relieve  the  excessive  training 
requirement  and  reduce  the  computational  complexity  of  the 
ACE  detector,  when  the  spatial-temporal  dimension  is  large. 
The  test  statistic  of  the  SI-PRao  test  has  been  asymptotically 
verified  to  be  independent  of  the  unknown  parameters  under 
the  null  hypothesis,  which  asymptotically  achieves  the  con¬ 
stant  false  alarm  rate.  Numerical  results  verify  that  the  SI- 
PRao  test  shows  better  performance  than  the  NPRao  test  and 
the  ACE  when  the  training  signals  are  limited. 

Appendix 

In  order  to  facilitate  the  derivation  of  the  SI-PRao  test,  we 
define  the  following  notations 


•  Gr  =  [au,ai,]T  =  [3?{a}  ,  3  {a}]T  denotes  the  sig¬ 
nal  parameter  vector,  where  ft  and  Sr  denote  the  real 
and  imaginary  parts,  respectively.  The  joint  probability 
density  function  (pdf)  under  Hq  and  the  pdf  under  Hi 
differ  only  in  the  value  of  Gr,  where  Qro  =  [0,0] T  and 
G  n  =  la.R,ai]T\ 

•  Gs  =  [A,  q^,  qf,  a#, &J ]  denotes  the  nuisance  pa¬ 
rameter  vector  with  =  vec  (3?  { AH  })  , &J  = 
vec  (Sr  {  A'7  }),  i[J_,  contains  the  diagonal  elements  in  Q 
and  the  real  part  of  the  elements  below  the  diagonal,  while 
q J  contains  the  imaginary  part  of  the  elements  below  the 

diagonal; 

r  -|  t 

.  G  =  Gj,,GTs  contains  all  unknown  parameters; 

r  t  ■'T  i T 

•  G  =  Gro ,  GSf)  denotes  the  ML  estimate  of  G  under 

H0. 

The  Rao  test  is  a  general  solution  to  a  class  of  parameter 
testing  problems.  It  is  often  simpler  than  the  GLRT,  and  is  also 
asymptotically  equivalent  to  the  latter.  A  detailed  discussion 
on  the  attributes  of  a  generic  Rao  test  can  be  found  in  [27], 
The  general  Rao  test  can  be  expressed  as  [27] 


din  f(G)  |T  r  ,  /-XI  din  f(G)  , 

OG,  ls=e  I  \  J\or,er  OG,  le=0’ 


(42) 


where 


C1"1  (*)]«„,«„ 

=  (l*r,er  (0)  -  I er,e.  (8)  l£0.  (8)  I e.,er  ( 8 ))  ,  (43) 

which  is  related  to  the  Fisher  information  matrix  (FIM)  [27] 


m  = 


ler,er  (8)  leT,eB  (8) 
U.,er{8)  I  e.fi.iP)- 


Under  both  hypothesis,  the  joint  pdf  of  the  test  signal  and 
training  signal  can  be  written  as 


f{8)  = 


A  1?+T 


exp  {  —  tr  ( Q  1T(A,  A))} 


(K+1)(N-P) 


where 


T(A,A)=- 


N-l  K  N-l 

j  E  £o  (n)  so11  (n)  +  Z  E  £k  (n)  ef  (n) 

n—P  k— 1 n—P 

:  (K  +  1)  (N  -  P )  ' 


£0  (n)  =x0  (n)  -  cts  (n) 


ek  (n)  =xfc  (n)  +  ^  aH  (P)xfc  (n~p)  ■ 


p 

x0  (n)  =x0  (n)  +  aH  (f)x o  {n  -  p) ,  (49) 

P=  1 
P 

1  (n)  =s  (n)  +  ^2  A H  ( V )  s  (n  —  p). 
p=  i 


(50) 
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(45)  gives  the  pdf  under  iTo  by  setting  a  =  0  and  the  pdf 
under  Hi  for  a  ^  0. 

From  (42),  the  derivation  of  the  SI-PRao  test  is  a  two- 
step  process.  The  first  one  is  to  obtain  the  ML  estimates  of 
the  nuisance  parameters  under  Ho,  and  the  second  one  is  to 
evaluate  the  related  terms  in  (42). 


A.  ML  Estimation  Under  Hq 

From  (45),  the  joint  pdf  under  H o  is  f(9)  with  a  =  0.  By 
taking  the  derivative  of  the  log  likelihood  In  f{9)  with  a  =  0 
with  respect  to  (w.r.t.)  Q  and  equating  it  to  zero  results  in  the 
ML  estimate  of  Q  as 

Qml  (A,  A)  =  T  (A,  A),  (51) 

where  T(A,  A)  is  given  in  (46)  with  a  =  0.  Substituting 
Qml  (A,  A)  into  In  f[0)  under  Hq  yields 

(K+1)(N-P) 

,  (52) 

The  ML  estimate  of  A  is  obtained  by  minimizing  |T  (A,  A) |. 
Note  that  ( K  +  1)  ( N  —  P)  T  (A,  A)  can  be  rewritten  as 

(K+  1)  (N-P)  T(A,A) 

=  (A)  +  R£  (A)  A  +  AHIlyX  (A)  +  AhR yv  (A)  A. 

It  can  be  shown  that  [12] 

T(A,A)  >  T(A,  A)  |a=a  ,  (53) 

where 

Aff(A)  =  -Rf,(A)R-1(A)  (54) 

and  RXX(A),  RW(A)  and  Rya:(A)  are  given  by  (31),  (32)  and 
(33),  respectively. 

Substituting  AH(A)  into  the  log  likelihood  In  f{6)  and 
ignoring  the  terms  independent  of  A  yield 

-  In  /  (a,  A,  Qml)  oc  ^  In  A  +  In  T  (a,  A)  ,  (55) 

where 

T  (A,  a)  |  oc  \RXX  (A)  -  Rfx  (A)  R"1  (A)  Ryx  (A)  |  ,  (56) 

and  the  symbol  oc  means  “proportional  to”. 

To  derive  the  ML  estimate  of  A,  we  need  to  evaluate 
T  ^A,  A^  .  Note  that  the  term  at  the  right  hand  side  of  (56) 
is  the  Schur  complement  of  [30] 


/  (\  A.  Qml) 


(e7r) 

A^tt  |t(A,A)| 
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From  (57),  (31),  (32)  and  (33),  it  is  shown  that  R(A)  can  be 
separated  into  A-dependent  and  A-independent  parts  as 

R  (A)  =  A_1R0  +  Rx,  (59) 

where  Ro  and  Rx  are  given  by  (22)  and  (23),  respectively. 
Note  that  Ro  and  R  r  are  positive  definite  with  probability  1 
if  N-P>  J(P+  1). 

To  differentiate  the  determinants  of  R  (A)  and  R,/:,;  (A)  w.r.t. 
A,  we  have  to  resort  to  the  following  Lemma. 

Lemma  1 :  Let  E  and  F  be  M  x  M  positive  definite  matrices. 
Then 


Id"  A- 


where  {£i}4  are  the  eigenvalues  of  E  ^FE  2. 
Proof: 


E  +  F  = 


M 

\  ' 


-1 


i=i  A(l  +  A^)’ 


d 

OX 


In  A_1E  - 


fi4‘" 

4>» 

4- 


E 


K 


A_1I 

M 


n  4 1 + &) 


.i= 1 


M 


=E 

i-1 

M 


d_ 

dX 


In 


1  +  A£i 

A 


-1 


a(i  +  a  eo’ 


where  we  have  used  the  fact  that  {A_1+£j}E  are  the 
eigenvectors  of  the  matrix  A_1I  +  E~iFE~5.  ■ 

By  applying  Lemma  1,  we  get 


(60) 


A  1I  +  E  2FE  2  )E2 
+  E~iFE”i 


(61) 


d_ 

dX 


In 


R(A) 


-1 


J(p+  1) 

E  A  (1  +  A  Hi) 


d  1 
~dX  ln 


R yy  iX) 


JP 


-1 


~iX(1  +  Xvi) ' 


(62) 

(63) 


where  /x,  and  tj,  are  eigenvalues  of  the  matrices  R0  2  R./\-Rn  2 
and  R0  2  Rx:j,R0  2 ,  respectively,  and  Ro,i/  and  R,A',y  are 
given  by  (24)  and  (25).  Following  that,  the  ML  estimate  of 
A,  denoted  as  A,  is  shown  to  be  the  root  of  (21).  As  a  result, 
the  ML  estimates  of  Q  and  A  are  given  by  (51)  and  (54), 
which  are  obtained  by  replacing  A  with  A. 


B.  Derivation  of  the  SI-PRao  Test 

From  (45),  the  elements  of  the  first  partial  derivative  of  the 
log  likelihood  In  /  with  respect  to  0,  in  (42)  are 


dhifjB) 

d0r 


dln/(9)  ‘ 

d\nf\e)  > 
dot.  i 


(64) 


with 


N- 1 


d  In  /  (0)  1  ^2  [§H  (n)  Q  1eQ  (n)  +  (n)  Q  !§  (n)] , 


daR 


N-l 


dlnf(e)  =j  e  j-gfl-  (n)  q  teo  (n)  _  £h  (n)  q  ig(n)]. 


dai 


^=p 


Since  the  term  [I  1  (0)]  g  g  is  related  to  the  FIM,  we  first 
find  relevant  entries  of  the  FIM: 


E 

E 

E 

E 


d2  In  /  (0) 

da\ 

d2 In  /  (0) 
da2 

d2 In  /  (0) 
daedal 

92  In  /  (0) 


xv-i 


=  -  t  E  gH  (n)Q  (4 , 


n—P 

N-l 


=  -  t  E  §ff  (n)Q  (n) , 


=0, 


=0, 


daida.R 
ler,es  (0)  =0, 

U.,0r  (0)  =0. 

As  a  result,  (43)  reduces  to 
[I-1  (0)] 


A 


er.er  N-i 

2  (n)Q^1s  (n) 

n—P 


1  0 
0  1 


Since 

91n/(0) 


\e=e 


daR 
N-l 

=  T  E  (n)  (n)  +  ^  (n)  Q  xs  (n) 

X  n=P  L 

d\nf(G)  , 

dai  '0=e 
.N-l 

~  3  ^2  (n)  Q_1x0  (n)  -  (n)  Q  xs  ( n ) 


A 


i=p 


(65) 


(66) 


(67) 


the  general  Rao  test  in  (42)  reduces  to  (20). 
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